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ABSTRACT 

Observations of 7-ray burst afterglows suggest that the correlation length of magnetic field fluctuations in 
the downstream of relativistic non-magnetized collisionless shocks grows with the distance from the shock 
to scales much larger than the plasma skin depth. We argue that this indicates that the plasma properties 
are described by a self-similar solution, and derive constraints on the scaling properties of the solution. For 
example, we find that the scaling of the characteristic magnetic field amplitude with distance from the shock is 
B oc D SB with -1 < sg < 0, that the spectrum of accelerated particles is dn/dE oc E~ 2 " SB+1 \ and that the scaling 
of the magnetic correlation function is (B,(x)B J (x+ Ax)) oc Ax 2s " (for Ax ^> D). We show that the plasma 
may be approximately described as a combination of two self-similar components: a kinetic component of 
energetic particles, and an MHD-like component representing "thermal" particles. We argue that the thermal 
component may be considered as an infinitely conducting fluid, in which case sb = and the scalings are 
completely determined (e.g. dn/dE oc E~ 2 and B oc D°, with possible logarithmic corrections). Similar claims 
can be made regarding non-relativistic shocks such as in supernova remnants, if the upstream magnetic field 
can be neglected. If self-similarity holds, it has important implications for any model of particle acceleration 
and/or field generation. For example, we show that the velocity-angle diffusion coefficient in diffusive shock 
acceleration models must satisfy D M/1 (p,D) = D~ l D^^{p/D) (where p is the particle momentum), and that a 
previously suggested model for the generation of a large scale magnetic field through hierarchical merger of 
current-filaments should be generalized. The generalization leads to qualitative changes in the model's major 
predictions, e.g. a constant merger velocity is predicted rather than increasing velocity approaching the speed of 
light. Finally, we point out that the self-similarity assumptions may be tested through the implied evolution of 
homogenous (time-dependent) plasmas, which may be accessible to direct numerical simulations: We predict 
that the inclusion (at the initial conditions) of a power-law spectrum of high energy particles, dn/dE oc E 2+l ', 
would lead to magnetic field evolution following B oc f~('i> +4 V 2 ('p+ 3 \ 
Subject headings: acceleration of particles — shock waves — gamma rays: bursts ■ 
mas — supernova remnants 
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1. INTRODUCTION 

Due to the low densities characteristic of a wide range of 
astrophysical environments, shocks observed in many astro- 
physical systems are collisionless, i.e. mediated by collective 
plasma instabilities rather than by particle-particle collisions. 
For example, collisionles s shocks play an importan t role in 
supernova remnants (e.g. Blandford & Eichl enl 19871). jets o f 
radio galaxies (e.g. Begelman e^JJ|1^9jyM^r^schi|j2003), 
gamma-ray bursts (GRB's, e.g. Zhang & Meszaros 2004), and 
the formation of the larg e scale structure of the Universe (e.g. 
lLoeb & WaxmanlbOOOh . Although collisionless shocks have 
been studied for several decades, theoretically and experimen- 
tally, in space and in the laboratory, a self-consistent theory 
of collisionless shocks based o n first princ iples has not yet 
emerged (see, e.g., comments in Krall 1997). 

Observations of GRB "afterglows," the delayed low en- 
ergy emission following the prompt 7-ray emission, provide 
a unique probe of the physics of collisionless shocks. Cur- 
rent understanding suggests that the afterglow radiation ob- 
served is the synchrotron emission of energetic, non-thermal 
electrons in the downstream of a strong collisionless shock 
driven into the surrounding interstellar medium (ISM) or stel- 
lar wind. These collisionless shocks start out highly relativis- 
tic, with shock Lorentz factor j s ~ 100 on time scale of min- 
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utes after the GRB, and gradually decelerate to 7., ~ 10 on 
a day time scale and 7, ~ 1 on a month time scale. This 
allows one to probe the physics of the shocks over a wide 
range of Lorentz factors. Afterglow shocks are highly "non- 
magnetized:" The ratio of magnetic field to kinetic energy flux 
ahead of the shock is very small, UB\/n\m p c 2 ~ 10~ 9 , where 
Ubi and ri\ are the magnetic energy density and particle num- 
ber density in the upstream rest-frame. This strongly suggests 
that the shock structure is determined by the upstream density 
and the shock Lorentz factor alone (e.g. Gruzinov 2001;j). We 
therefore adopt the assumption that the shock structure ap- 
proaches a well defined limit as UB\/n\m p c 2 — ► 0, and that 
GRB afterglow shocks are approximately described by this 
limiting solution (see § 12. ll for a detailed discussion). 

It may be noted here that the upstream density can 
be eliminated from the problem by measuring time in 
units of the (shock-frame proton) plasma time, ui' 1 = 

(4-nj s nie 2 /m p 'y s y~ 1 ' 2 = (A-Kn\e 2 /m p )~ l l 2 , and measuring dis- 
tances in units of the corresponding skin depth, l sc / = c/ui p . 
The shock is then completely specified by the dimensionless 
parameter 7 s v s /c (and the dimensionless mass ratio m e /m p ; 
the upstream pressure is assumed negligible). In this sense, 
GRB shocks may be considered "simple." 

The synchrotron model of GRB afterglows requires a strong 
magnetic field and a population of energetic electrons to 
be present in the downs tream. Optical observations (e.g. 
Zhang & Meszaros 2004), the clustering of explosion ener- 
gies {Frail et al. 1 1200 1|) . and the observed X-ray luminosity 
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jFreedman & Waxmanl200illBerger et a l. 2003) suggest that 
the fraction of post-shock thermal energy density carried by 
non-thermal electrons, e e , is large, e e « 0.1. The fraction of 
post-shock thermal energy carried by the magnetic field, eg, 
is less well constrained by observations. However, in cases 
where eg can be reliably constrained by multi waveband spec- 
tra, values close to equipartition, en ~ -01 to 0.1, are inferred 
(e.g. lFrail. Waxman & Kulkarnil2000l) 3 . 

Near equipartition magnetic fields may conceivably be pro- 
duced in the collisionless shock driven by the GRB ex- 
plos ion by elec tromagne tic (e.g. Weibel-like) instabilities 
(e.g.lBlandford & Eichle rlll987t lGruzinov & Waxman _lfT999t 
Medv edev & Loebll999HWiersma & Achterbergl2004l) 4 . The 
main challenge associated with the downstream magnetic 
field is related to the fact that in order to account for the 
observed radiation as synchrotron emission from accelerated 
electrons, the field amplitude must remain close to equipar- 
tition deep into the downstream, over distances ~ lQ 10 l sc f. 
At t ~ 1 d the magnetic field must be strong through- 
out the (proper) width A ~ 2j s ct ~ 10 17 cm while l sl j ~ 
lO^ni/lcm" 3 )" 1 / 2 cm. This is a challenge since electromag- 
netic instabilities are believed to generate (near-equipartition) 
magnetic fields with coherence length L ~ l s j, and a field 
varying on such scale is expected to dec ay within a few 
skin-depths downstream (Gruzinov 2001a). This suggests 
that the correlation length of the magnetic field far down- 
stream must be much larger than the skin depth, L ^> l st /, 
perhaps even of the order of the distance f rom the shock 
(Gruzinov & Waxman 1999; Gruzinov 2001a). 

Growth of the characteristic length scale by many orders 
of magnitude, from L ~ l s j to L ^> l s d, is a strong indication 
of self-similarity. In regions where L 3> l s j, which also im- 
ply that L is much larger than the Larmor radius of thermal 

protons (L 3> Ri.rh ~ 7s'«pC 2 /eB ~ hd^s ), it is reasonable 
to assume that L is the only relevant length scale, and self- 
similarity is expected (see §|3]for a detailed discussion). The 
main goal of the current paper is to introduce and formulate 
the assumption of a self-similar collisionless shock structure, 
and to study some of its consequences. 

Although our analysis is motivated by GRB afterglow ob- 
servations, it may be relevant also for non-relativistic colli- 
sionless shocks, such as shocks in young supernova remnants 
(SNRs). In the past few years, high resolution X-ray obser- 
vations have provided indirect evidence for the presence of 
strong magnetic fields, > 100 /iG, in the non-rel ativistic (v., ~ 
few x 1000 km s" 1 ) shocks of young SN Rs (see Bamb aet al.l 
l200l IVink & LamindEOQl IVolk et alJl200l . These fields 
extend to distances D > 10 17 cm ~ 10 lo hd downstream, and 
possibly even > 10 16 cm upstream, of the shock. In resem- 
blance to GRB's, such strong magnetic fields cannot result 
from the shock compression of a typical interstellar medium 
(ISM) magnetic field, B\ ~ few /jG. In SNRs, the discrep- 
ancy is somewhat less severe, UB\/n\m p v 2 s < 10~ 4 , and the 
possibility that these magnetic fields are related to the large 

3 Eichler & Waxman 1 2005 ) have pointed out that observations determine 
e e and es only up to a factor /, the fraction of electrons accelerated, where 
m e /m n < f < 1. Howeve r, it is expected that / is not very small, f > 1/10 
lEichler & Waxmarl2005l) . 

' It should be noted here that if a hydrodynamic description of the shock 
were applicable, the shock would ha ve been stable and no m agnetic fields 
could have been generated 1 Gruzinov 2000; Wane et al. 2002). A possible 
caveat is the presence of large inhomogeneities in the upstream plasma (ISM), 
which may drive turbulence in the downstream plasma, which may generate 
strong magnetic fields with large correlation length 



scale ISM fields cannot be ruled out. If the ISM magnetic 
fields can be neglected, this suggests that these shocks too 
may have a self-similar nature. Henceforth, when discussing 
non-relativistic shocks, we assume that this is indeed the case. 

The non-thermal energetic electron (and proton) population 
is believed to be produced by the diffusive (Fer mi) shock ac- 
celer ation (DSA) mechanis m (for reviews se e lDrurvlll983l 
iBlandford & Eichlerifl987t |Malkov & Drurvll2001l) . Accel- 
eration of charged particles to high, non-thermal energies 
is a ubiquitous phenomenon in both relativistic and non- 
relativistic collisionless shocks. The accelerated particles are 
estimated to carry a considerable part of the energy: elec- 
trons alone carry ~ 10% o f the thermal energy in GRB exter- 
nal shocks (Eichler & Waxman 2005) and ~ 5% of the ther- 
mal energy in SNR shocks ( Keshet et al. 2004, and references 
therein); and at least 10% of the energy in SNR shocks must 
be converted into relativistic proto ns if these shocks are re- 
sponsible for Galactic cosmic rays (iDrurv et alJl!989l) . This 
has several important implications. The accelerated parti- 
cles are likely to have an important role in generating and 
maintaining the inferred magnetic fields. This conclusion 
is supported also by the evidence of strong amplification of 
the magnetic field in the upstream of GRB afterglow shocks 
(Li & Waxman 2006), which is most likely due to the stream- 
ing of high energy particles ahead of the shock. Since the 
high energy particles are likely to play an important role in 
the generation of the fields, a theory of collisionless shocks 
must provide a self-consistent description of particle acceler- 
ation, which depends on the scattering of these particles by 
magnetic fields, and field generation, which is likely driven 
by the accelerated particles. 

The search for a self-consistent theory of collisionless 
shocks has led to extensive numerical studies. Particle in 
cell ( PIC) simulations were performed in one dimension 
(e.g. iDi eckmann et all l200 6|). in two d imensions (2D; e.g. 
Wallace & Epperleinl ll991t iKatol 120051) . in two spatial and 
three velocity d imensions (2D3V; e.g. iGruzinovl 1200 lab: 
iMedvedev et afll2005l and in the past few years also in t hree 
dimensions (3D; e g.lSilva et all2003tlNishikawa et al.l2003l: 
[ Frederiksen et al.l 120041: Uaroschek et alJ 120041: iSpitkovskvl 
2005). Such simulations have provided compelling evi- 
dence that transverse, electromagnetic (Weibel-like) insta- 
bilities generate near-equipartition magnetic fields in pair 
(e + e~) plasma, and eg > m e /m p magnetic fields in ion-electron 
plasma. However, 3D simulations are limited to very small 
simulation boxes, and at present can reliably probe small 
length scales no larger than ~ 100 electron skin-depths, and 
short time scales no longer than ~ 100 electron plasma times. 
Hence, there is only preliminary evidence for the existence of 
collisionless shocks in 3D, and only in a pair-plasma (Simula- 
tions of ion-electron plasma are forced to employ an effective, 
small proton to electron mass ratio, m p /m e < 20 with present 
computational resources, and the preliminary results thus ob- 
tained are not easily extrapolated to more realistic mass ra- 
tios). Obviously, the question of field survival and correla- 
tion length evolution on length scales 3> l s d are not yet an- 
swered. Similarly, highly energetic particles cannot be con- 
tained in the small simulation boxes used, so Fermi-like accel- 
eration processes are suppressed. It is important to note, that 
some published results are based on PIC simulations in stages 
where the boundary conditions strongly modify the plasma 
evolution. For example, claims that the magnetic fields decay 
slowly or saturate at some finite level remain questionable, 
until verified by simulations with sufficiently large simulation 
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boxes. A discussion of 3D PIC simulations and their physical 
implication appears in Appendix $ lAl 

In §|2] we lay the basis for our analysis of aftergl ow ( rel- 
ativistic) and SNR (non-relativistic) shocks. In § 12.11 we 
discuss our assumption that these shocks are highly "non- 
magnetized," i.e. that the shock structure approaches a well 
defined limit as Ubi /n\tn p v 2 — > and that the shocks obse rved 
are approximately described by this limiting solution. In § 12.21 
we present the governing equations and discuss their depen- 
dence on dimensional parameters, which is relevant for the 
discussion of self-similarity. In particular, we demonstrate 
that when distances are measured in units of l s d, th e sho ck 
structure depends only on 7,v,/c (and m e /m p ). In § I2.3l we 
clarify the notion of "shock structure:" Since the electromag- 
netic fields and particle distributions fluctuate with time (at 
any given point downstream, and perhaps also upstream of 
the shock), the "stationary shock structure" is given by the 
correlation functions of the fluctuating quantities (which are 
expected to depend only on the distance from the shock). We 
define a notation that is useful for the discussion of the self- 
similarity assumption. 

In §|5]we introduce and formulate the assumption of a self- 
similar structure in the downstream of non-magnetized colli- 
sionless shocks, derive several scaling relations of the physi- 
cal quantities, and discuss some of the physical implications. 
One of the conclusions of §|3]is that the plasma may be ap- 
proximately described as a combination of two self-similar 
components: a kinetic component of energetic particles, and 
an MHD-like component representing the bulk, "thermal" par- 
ticles. The MHD-like component is discussed in § 0] We 
argue that this component may be treated as an infinitely con- 
ducting fluid, and show that this leads to a complete determi- 
nation of the scaling laws. 

In §|5]we present various extensions of the analysis. Self- 
similarity is studied in the upstream of non-magnetized col- 
lisionless shocks and in homogenous time-dependent plas- 
mas, which may be more accessible to simulations than (non- 
homogeneous) collisionless shocks. In §[6] we discuss some 
of the implications of the self-similarity assumption to mod- 
els of diffusive particle acceleratio n, and to the phe nomeno- 
logical model, suggested by Medvedev et al. (2005), of field 
generation through hierarchical merger of electric current fil- 
aments. The latter model is generalized, and shown to follow 
the self-similar scalings. Our main results and conclusions are 
summarized in §0 

2. NON-MAGNETIZED COLLISIONLESS SHOCKS 
2.1. Do non-magnetized collisionless shocks exist? 

In general, the parameters determining the structure of a 
collisionless shock (for a given homogenous upstream com- 
position which we assume to be an electron-proton plasma) 
are the shock (four) velocity 7 s v s , and the upstream density 
ri\, pressure P\ and magnetic field Bi. We can formally ask 
what the limiting configuration is when Bi and Pi approach 
zero. We assume that in this limit there is some shock so- 
lution, to which we refer as non-magnetized and strong. It 
is important to stress that by "non-magnetized" ("strong") we 
refer to the null effect of the value of the upstream magnetic 
field (pressure) on the stationary configuration of the shock, 
regardless of the dynamical role of these parameters in the 
generation of the shock. 

We next ask how small should B; and Pi be so that the 
solution is no longer affected by their value. The dimension- 



less parameters ub = UBi/nim p v 2 ,ap = P\/n\m p \ 2 measure 
the relative contribution of the magnetic field and thermal en- 
ergy to the energy flux in the shock frame (which is constant 
and thus equal to its value in the far upstream), 

T 0: ~ jjn i m p v 2 ( l + a B + <Jp)v s , ( 1 ) 

where we assumed that the shock normal (the direction of 
the flow) is in the z direction. It is reasonable to assume 
that if as and o> are much smaller than unity, the magnetic 
field and pressure do not affect the shock structure. Since 
o~b ~ (lsd I Ri,th) 2 , the condition erg <C 1 implies that the ther- 
mal protons in the shock frame are not affected by the com- 
pressed upstream magnetic field on the dynamical distance 
~ l s d (the cyclotron time ~ ^ s m p c/eB is longer than the dy- 
namical time ~ u!~ l as long as o~b < (c/v,) 2 ). 

Simulations of relativistic shocks with various values of 
the upstream magnetic field (for a shock with Lorentz factor 
T = 30 in electron-positron plasma) are reported in Spit kovskvl 
(2005), and a typical value of o~b ~ 10~ 2 is found to distinguish 
between the regimes of magnetized and non-magnetized 
shocks. Note, that this result may change considerably for 
electron-ion plasmas or if a population of accelerated parti- 
cles is included, and may be affected by the transverse (per- 
pendicular to the flow) boundary conditions (see §|A). There 
is some debate regarding our assumption that the magnetic 
field i n the upstream of GRB external shocks is insignificant. 
Lvubarskv & Ei chlerl ( 120051) claim that the initial stages of 
the generation of the shock depend on the upstream magnetic 
field, because it can isotropize the protons faster than the mag- 
netic fields generated in the shock can. Note, however, that 
even if the upstream magnetic field does play an important 
dynamical role in the evolution of the shock, this does not im- 
ply that the steady-state solution depends on the exact value 
of the upstream ma gnetic field. 

Simil arly, iBelil (120041 120051) and iMilosavlievic & Nakarl 
(2005b) claim that the cosmic ray precursor of the shock can 
amplify the upstream magnetic field significantly. Assuming 
that the magnetic field is amplified by many orders of mag- 
nitude, it is reasonable to assume that its initial value is not 
important. 

2.2. Governing equations and dimensional considerations 

Henceforth, we assume that strong, non-magnetized colli- 
sionless shocks do indeed exist. We consider a quasi steady- 
state planar shock. First, we consider the equations that gov- 
ern the full plasma distribution function, / Q (x,p,f), and are 
valid globally, i.e. both upstream and downstream. Subscripts 
a e {e,i} denote electrons and ions, respectively. Approx- 
imate equations, which are valid only in the far downstr eam 
and f or wh ich self-similarity is assumed, are discussed in § 13.31 
and § |4~T1 

The flow is governed by Vlasov's equation, 

d,f Q + v(p) ■ Vf a + q a (E+ I v x • V p / Q = 0, (2) 
and Maxwell's equations, 

4vr 1 
VxB = — j+-d,E, 

c c 

VxE = --d,B, (3) 

c 

where the electric current is related to f a through 

j = X^ Q M 3 PaV(p)/a(p). (4) 
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The velocity of the particles is given in terms of the momen- 
tum by 

v(p) = c^=i_. (5) 

As usual, the equations V • B = and V • E = 4np are assumed 
to hold at some initial time, and are therefor preserved by the 
other equations at all times. 

The shock is completely defined by the far upstream bound- 
ary conditions. Written in the shock frame, these conditions 
are 



E(D -> -oo) = B(D -> -oo) = 0, 
f a (D — > -oo) = 7. ! «ic$ 3 (p-7. 5 v J TO Q z), 



(6) 



where D is the distance from the shock front in the direction 
of the downstream (D is negative in the upstream). Note that 
the equations are independent of the frame of reference. 

In order to highlight the dimensional dependencies of the 
Vlasov-Maxwell equations, we express the physical quanti- 
ties A = t ,x . . . in term of dimensionless variables A, 

Isd - , - 
t= — t, X = l s dX, p= p, h p, 

v., m„ 



and 



j = ev. s 7 J n 1 j, 
U = 



(7) 



3 J a i 



(8) 



where 



\-Kn\e 1 47r«2e 2 



m„ 



JsTTlp 



C 

kd = , 

Pth = m p j s v s 



(9) 



are the characteristic ion plasma frequency (squared), ion 
skin-depth and characteristic momentum of thermal particles 
in the downstream. Inserting these into Eqs. and Q, we 
obtain 



V s m a /m p 



E + -xB )-V p f a =0, 



VxB = j + ^a,-E, 
V x E = -d f B 



and 



with 



J = 



d P— fa, 

v., 



(10) 



(11) 



V 

Vv 



l i v ; -? 



The boundary conditions at upstream infinity (written in the 
shock frame) can similarly be written in dimensionless form, 

E0 — > -oo) = B(D — > -oo) = (13) 

and 

/ Q 0^-oo) = (S 3 (p-z), 

(14) 

where D = D/l s ^. 

Measuring distances in units of l s j, we have thus arrived at a 
set of dimensionless equations, Eqs. dl0t -( ll4t . which depend 
only on 7 s v s /c (and on m e /m p ). 

2.3. Stationary shock structure 

The electromagnetic fields and particle distributions at any 
given point behind (downstream of) the shock fluctuate with 
time. For each set of values of the shock parameters, n\ and 
7sV s , there is a large ensemble of time-dependent "specific 
solutions," with specific temporal and spatial dependence of 
particle distribution functions and electromagnetic fields. We 
assume that the averages and correlation functions of the fluc- 
tuating quantities depend only on the distance from the shock, 
and are identical for all specific solutions (in the limit that the 
size of the shock plane is infinite). 

Consider for example the particle distribution function 
/ Q (x = Dz+x^,p,f) or magnetic field B(x = Dz+x^,f), where 
we have separated the x dependence to dependence on the 
distance from the shock front, using z as the direction of the 
shock normal, and onxi, a two dimensional vector perpen- 
dicular to z. While f a depends on t and x^, we assume that an 
average of f a over planes perpendicular to the shock normal 
is independent of t and x^: 

(/a)(p;£>) = lim ~K I d 2 x' ± f a (Dz+x ± +x^,p,f) . 

r - , °° 7rr* J x ' ±<r 

(15) 

Similarly, we assume that correlation functions of /, B and 
E may be defined, which are independent of t and xj_, e.g. 

B, 7 (Ax,Af,D) = 

1 



lim 



, , d 2 x' ± 



(12) 



Bj [Dz+x_L+x' ± ,t] Bj [Dz + xj_ + Ax+x^_,f + Af] . 

(16) 

The stationary shock structure is described by (f a ) and 
by the (infinite) set of correlation functions. Note, that the 
Maxwell- Vlasov equations may be converted to an (infinite) 
hierarchy of equations for the correlation functions. 

It is useful for the analysis that follows to define the two- 
dimensional power spectrum of the magnetic field in a plane 
parallel to the shock front. We define the 2D power spectrum 
at distance D as 

B ( p(k ± ,D) = (27T)" 1 J d 2 x' ± B,/xl,zV^< , (17) 

where B u {x'^ , D) = B u { Ax = x' x ,At = 0,D). The plane av- 
eraged magnetic energy density is related to the correlation 
function by 

(%)(D) = (87r)- 1 ^(27r)- 1 f d 2 k ± Bf°\k ± ,D) 

= (87r)- 1 ^B // (Ax = 0,Af = 0,D). (18) 
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Note that the plane-averaged magnetic field vanishes. As 
a consequence of the planar symmetry, a non-zero averaged 
magnetic field would be oriented along the z-axis. As its di- 
vergence vanishes, the magnitude of (B) would be constant 
and thus equal to its far upstream value which is zero. The 
root mean square (rms) magnetic field, 



B(D) = 



^fi,,(Ax = 0,Ar = 0,D) 



1/2 



(19) 



serves as a measure of the characteristic magnetic field ampli- 
tude at a distance D from the shock (note, (Ub) = B 2 /8ir). 

3. SELF-SIMILARITY IN THE DOWNSTREAM 

As discussed in the introduction, afterglow observations 
suggest that the characteristic length scale L for variations in 
the magnetic field becomes much larger than l s d, L^S> l s d, at 
distances D ^> /„/ downstream of the shock. We assume here 
that in the limit L/I s a — > oo, L becomes the only length scale 
relevant for the evolution of the plasma, which implies self- 
similarity. There is no proof that self-similarity will be present 
whenever the characteristic length scale diverges (or becomes 
infinitesimal). However, the self-similarity assumption is 
known to be valid for many physical systems in which s uch 
divergence occurs [see, e.g., Zeld ovic h & Raiz erl (U96 8) for 
self-similarity in hydrodynamics, and Kadanoff et al. (1967) 
for self-similarity in critical phenomena]. In order to clarify 
the reasoning behind thi s ass umption and its implications, we 
first briefly describe in § !3.1l an example of a physical system, 
where the presence of a diverging characteristic leng th sc ale 
leads to self-similar behavior. We then formulate in § l3.2l the 
self-similarity assumption for the downstream of collisionless 
shocks. The s epara tion of the plasma into two compone nts is 
discussed in § 13.31 The scaling relations are deri ved in § 13.41 
and some of their implications are discussed in § 13.51 

A note is in place here regarding non-relativistic collision- 
less shocks. As mentioned in the introduction, self-similarity 
may be expected in the non-relativistic shocks observed in 
young SNR's. The analysis is similar for relativistic and non- 
relativistic shocks, and in particular leads to the same scaling 
relations and implications. The main difference is that in the 
non-relativistic case there is a physically relevant length scale, 
which is different than the skin depth l sc i'- The Larmor radius 
of mildly relativistic protons, Ri ~ m p c 2 /eB ~ R-L,th(Vs/c)~ l ~ 
(v s /c)~ l lsd, where Ri^h is the characteristic Larmor radius of 
thermal protons. In other words, there is a small parameter 
in the problem, (v s /c), which is likely to have a physical im- 
plication. Hence, although in this case too we expect self- 
similarity as L/l s d — > oo, we expect the self-similar solution 
to provide a good approximation for L ^> (y s /c)~ l l sc i rather 
than for L ^> l st j, and the scaling relations of the particle dis- 
tribution functions to be valid for momenta p > m p c. 

3.1. An example of self-similarity 

Consider a macroscopic system in thermodynamic equilib- 
rium, exhibiting a phase transition at some temperature T c . 
Near the critical temperature, T ~ T c , the correlation length 
L of fluctuations within the system diverges, L(5T) — > oo as 
ST = T - T c — ► 0. It is commonly assumed that in this limit L 
becomes the only "relevant" length scale, and that the micro- 
scopic length scales, e.g. the inter-particle distance d, become 
irrelevant. 

To clarify this assumption, let us consider the correlation 
function f(r) of some fluctuating quantity. / is a function 



of 5T and of the various parameters defining the thermody- 
namic system (e.g. d). Assuming L diverges monotonically 
with ST we may replace the dependence on ST with depen- 
dence on L, f(r;L,d,Ci) where c,- represent the parameters 
defining the thermodynamic system. We may assume that 
these parameters contain no parameter with the dimensions 
of length, since any such parameter q may be replaced by 
a dimensionless parameter, q jd. Using the dimensional pa- 
rameters, a constant A f with the same dimensions of / may 
be constructed, and / may be written as / = Aff where / is a 
dimensionless function. Since / is dimensionless, it must be 
a function of only dimensionless combinations of {r,L,d,Ci}. 
Thus, f =Aff(r/L;d/L,dj) where c,- represent the dimension- 
less combinations of c,. As L/d diverges it is assumed that d 
becomes "irrelevant," in the sense that / approaches a limit 

f(r;L,d,Ci)=A f f (^;-,Ci ) —^Ai I - 



d/L^O " \L / 

(20) 

Note, that it is not assumed that, as d/L — > 0, / approaches 
a limit independent of d/L, but rather that in this limit it 
has a scale-free (power-law) dependence on d/L. In the for- 
mer case, / depends on the dimensional parameters A f and 
L alone, and in the latter on A fd~ Sf and L. In both cases, the 
only length scale that may be extracted from / is L. Simi- 
lar arguments lead to the conclusion that L has a power-law 
dependence on ST, 



6T(L,d,a) = T c ST 



-iCi 



d/L^O 



hied- (21) 



In general, the self-similarity assumption may be stated as 
the assumption that all the properties of the system at some 
ST = STi, corresponding to L\ = L(ST\), are identical to those 
of the system at ST = 6T2, corresponding to Li = L{5T2), up 
to a scaling transformation. That is, any function describing 
some properties of the system, e.g. f[r, 8T(L)], is identical at 
L\ and L2 up to a scaling of r and /, 

'L2 



f(r,L 2 ) = 



f 



L2/L1 



(22) 



for some Sf. The function / is termed self-similar, as it is 
similar to itself at different L values. Comparing Eqs. J22I 
and < !20b it is clear that the assumption that d is "irrelevant" is 
equivalent to the assumption that / is self-similar. 

It is clear from the above example that the assumption of 
self-similarity provides powerful constraints on the proper- 
ties of the system, and that a complete characterization of its 
properties requires determination of the similarity exponents 
{s f ,s T }. 

3.2. The self-similarity assumption 

Consider the plasma at a distance D l s d downstream of 
the shock, where the field correlation length is assumed to 
satisfy L ^> l s d- Assuming L/l s d diverges with D, we expect 
the structure of the shock to become self-similar. That is, we 
expect the averaged distribution function (f a ), the rms mag- 
netic field B and the (infinite) set of correlation functions to 
be self-similar: 

B(L)=(J-) B(Lo), (23) 



(fa)(P,L) = 



(fa) 



(L/LoY- 



(24) 
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B, 7 (Ax,Af,L) = 



La 



2sb 



Ax 



At 



L/L Q : (L/L y 



-,Lq 



, (25) 



and so on, where Lq is a reference length scale. Note, that we 
have replaced the dependence on D with a dependence on L, 
assuming that L diverges monotonically with D. 

The values of the similarity exponents are not determined 
by the requirement of self-similarity alone, and must be de- 
rived based on additional arguments. In what follows, we de- 
rive constraints on the similarity exponents {s r , s p , s B ...} us- 
ing the Maxwell-Vlasov equations. Since self-similarity ap- 
plies to averaged quantities, e.g. (/„) and Bh, rather than to 
specific solutions, this derivation should be based on the equa- 
tions for the averaged quantities, which are obtained from the 
Maxwell-Vlasov equations. Such a derivation is, however, 
rather cumbersome. Identical constraints on the similarity ex- 
ponents may be derived in a simpler way by assuming that 
the specific solutions form a scalable family in the sense that 
if {/ o (i,p,0,B(x,0} is a solution of the Maxwell-Vlasov 
equations, then the scaled functions 



fa O^iPi = £ 1 fa ( ~Z j ~T7~ , ~rr 

u f" i* 



(26) 



also constitute a solution (at least approximately). It is 
straightforward to verify that the self-similarity requirements, 
Eqs. (I24M25> . are automatically satisfied under this assump- 
tion, by substituting £ = L/Lo in Eq. J26b and recalling that the 
averages and correlation functions are identical for all specific 
solutions. Under the assumption that the solutions are scal- 
able, we can use the Maxwell-Vlasov equations directly (in- 
stead of the equations for the correlation functions) in order 
to derive the constraints o n the similarity exponents. Before 
doing so, we address (in § 13 .3i the separation of the plasma 
into two components. 

It is important to emphasize here, that although the con- 
straints derived under the assumption that the solutions are 
scalable (in the above sense) are identical to those derived 
from the self-similarity requirements, Eqs. J24M 251. it is 
not obvious that the self-similarity requirements indeed im- 
ply that the specific solutions are scalable. 

3.3. Two plasma components 

As the plasma flows away from the shock, most of the parti- 
cles remain "thermal," 5 i.e. most protons carry a momentum 
~ p,h (and most electrons carry a momentum (m e /m p ) u p t h 
with v > 0). This implies the existence of a length scale 
Rh,th{L) = pthc/eB oc Lr SB , Rui,/L oc ZT (1+ * b) . Obviously, in or- 
der for L to be the only relevant scale, we must have Sb > — 1. 
Combined with the requirement that the field energy density 
does not diverge we have 



-1 <s B <0. 



(27) 



Our assumption that L is the only relevant length scale im- 
plies that Ri ,i, is irrelevant. However, Ruh is obviously rel- 
evant for the description of the microscopic motion of indi- 
vidual "thermal" particles. This apparent contradiction may 

5 The term "thermal" is used for convenience to describe particles with 
kinetic energy ~ (7, - l)m p c 2 . In what follows, it is not claimed or assumed 
that the energy distribution of these particles is thermal i.e. a (relativistic) 
Maxwellian distribution. 



be resolved only if the microscopic motion of "thermal" par- 
ticles, i.e. of particles with momenta p such that pc/eB <C L, 
is unimportant. These particles should therefore be described 
by effective equations, where the microscopic particle motion 
is unimportant. We leave the discussion of the "thermal" par- 
ticles, to which we refer henceforth as the "fluid" particles, to 
§0] and note here only several points which are important for 
the discussion that follows. 

First, the scaling of f a derived in the previous section does 
not apply to all particle momenta, but rather to the large p, 
p ^> Pthi^s/c)" 1 , behavior of f a . In particular, the Vlasov 
equation in this region, where v(p) w cp, can be written as, 

d,f a + cp ■ Vf a + q a (e+ iv x • Vp/ Q = 0. (28) 

Second, the total electric current is the sum of the electric 
currents carried by the high energy particles, j A , and by the 
"fluid," j F , j = j h +j F . Hereafter, h and F subscripts refer to 
the high energy and to the fluid particles, respectively. The 
electric current j h is given by 



J/, : 



/ 



d^cpfjp). 



(29) 



p>i p eBLjc 



Here, £ P (L) is a dimensionless function that determines the 
threshold momentum between accelerated and fluid compo- 
nents. It must satisfy £ p > (v s /c)~ l p t hc/eBL= (Vj/c) _1 /?£, jf ft/L 
in order to ensure that the accelerated component is not af- 
fected by the thermal scale. It must also tend to zero as 
L diverges, as self-similarity is expected for all scales S> 
(v s /c)-\R L ,th. 

Solutions of the type we have arrived at, where the self- 
similar solution describes the evolution in some part of (x, p, t ) 
phase space while the other part is described by a differ- 
ent solution, are usually ca l led "second type solutions" (e.g. 
IZeldovich&Raizerl 119681 IWaxman & Shvarts I 119931) . In 
second type solutions, the similarity exponents can not be 
determined by global conservation laws. For example, one 
may have argued that a relation between sj and s p may be de- 
rived by requiring the integral over momenta of (/) , given by 
Eq. (1241 . to be equal to the (downstream) particle density, 

n a = Jd 3 p(f a ){ V ,L)=(J-^ ' ' J d 3 y(f a )(y,L () ), (30) 

which would imply 3s p + s/ = 0. Such a constraint can not 
be derived, however, since the scaling relation, Eq. 124\ . used 
for obtaining the second equality of Eq. J30i does not hold 
for small values of y. In fact, as we show below, the func- 
tional form of / derived by the self-similarity arguments is 
such that the integral on the rhs of Eq. J30l > diverges at small 
>'. The number of particles described by the self-similar solu- 
tion does not diverge, however, since the integration extends 
only down to y = £, p (L)(L/Lq) 1 ~ Sp eBLo/c (note, that this con- 
strains the dependence of £ p on L). This is analogous to the 
divergence of energy in s econd-type self-similar s olutions of 
hydrodynamic flows (e.g. Waxman & Shvarts 1993). 

The particles with p < £ p (L)eBL/c are described by a so- 
lution different than the self-similar solution describing the 
higher momenta, "accelerated," particles. As in any second 
type self-similar solution, it is necessary that the non self- 
similar part of the solution does not affect the self-similar 
part (this requirement often determines the simi larity expo- 
nents of second type solutions, e.g. Zeldovich & Raizer 
119681 IWaxman & Shvarts lfl993l) . In our case, the particles 
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at p < £ p (L)eBL/c may affect the higher momentum particles 
only through their contribution to the electric current. This 
implies that the current contributed by p < ^ p (L)eBLjc parti- 
cles must either be negligible or scale with L in a similar way 
as the current of the higher momenta particles does. For the 
fluid current, this requirement implies that either jp <C jh or 
Jf(L)/ jh(L) oc L . In addition, this requirement implies that 
the integral on the rhs of Eq. ( I29l > converges (as otherwise the 
current would be dominated by p < ^ p {L)eBL/c particles). 

The following point should be emphasized here. We have 
implicitly assumed above that the high momenta, accelerated 
particles are dynamically important in the sense that their 
electric current makes a considerable contribution to the to- 
tal electric current, j h ~ j. However, self-similar solutions 
where the accelerated particles do not contribute to the cur- 
rent, jf, <C j, are possible in principle. In this case, the growth 
of the magnetic field fluctuation correlation length should be 
driven by the fluid, and the accelerated particles can be treated 
as "test-particles," which do not affect the flow. As mentioned 
in the introduction, such a picture is unlikely, due to the fact 
that the accelerated particles carry a significant fraction of the 
energy and due to the evidence that they play a role in gener- 
ating magnetic fields in the upstream plasma. 

3.4. Scaling relations 

Let us consider first the scaling of L with D. Since the char- 
acteristic length scale for changes in L is L, we must have 
LocD: Requiring the scale for changes in L to be proportional 
to L may be written as d(\ogL)/dD oc l/L which implies 
L oc D. It should be noted that this assumption is not equiv- 
alent to assuming that the only length scale relevant for the 
evolution is D. Under such an assumption, d(\og L) j dD = s/D 
implying L = D s , which allows the possibility L/D — > as 
D^oo. 

Next, we consider the scaling of time, i.e. the value of s t . 
Substituting the scaled solutions, Eqs. i26\ . into the Valsov 
equation, Eq. d28i . one finds that the various terms in the 
equation scale differently with £. In order for the scaling of 
the temporal and the spatial derivative terms to be similar, one 
must require s t = l. However, one can not conclude from this 
that s t = 1 is the only value allowed, since it is possible that 
one of the first two terms, d, f or v • V/, becomes negligible as 
L — > oo, and may be neglected altogether in the Vlasov equa- 
tion. Nevertheless, we conclude that s t = 1 is required based 
on the following arguments. The case s t < 1, implying that 
cT jh — > as L — > oo where T is the characteristic time for 
variations in the physical quantities, is ruled out since it would 
imply that at large distances from the shock the electric field 
is much larger than the magnetic field (E/B — > oo as L — > oo), 
which is inconsistent with the synchrotron model of afterglow 
emission. Assuming s, > 1, cT/L — ► oo asL — * oo, would im- 
ply that at large distances from the shock the magnetic fields 
are approximately static (velocities of the field lines of the or- 
der L/T) in the shock frame. Once the condition L/T <C v s 
would be reached, the thermal protons would gyrate around 
the static field lines (the gyration radius, R.L,th, is much smaller 
than L, and thus also the gyration time, Tl ,i, ~ Rl,iIi/vs, is 
much shorter than T), in contradiction with the fact that the 
bulk downstream fluid moves with velocity vj < v 5 (Note that 
7j = (1 - v^/c 2 ) -1 / 2 ~ 1). We conclude that we must have 
s f = l. 

The scaling of the momentum, i.e. the value of s p , is deter- 
mined by comparing the momentum derivative term with the 
spatial (or temporal) derivative term in the Vlasov equation. 



Substituting the scaled solutions, Eqs. (1261 . into the Vlasov 
equation, Eq. (I28> . and requiring all terms to scale similarly 
implies s p = \+sb- One implication of this scaling is that the 
Larmor radius scales as L (in fact, the entire trajectory of each 
particle scales with L). 

Finally, we derive the relation between Sf and sb- The cur- 
rent provided by the accelerated particles is 

= ^ +3ta+1) j A (|0. (3D 

We have not shown here explicitly the lo wer limit of the in- 
tegration over p, since, as discussed in § 13.31 the current in- 
tegral must converge at small p. Substituting the scaled cur- 
rent, and the scaled magnetic field, into Maxwell's equation, 
V x B = 47rc" I j + c" I <9,E, we find that the V x B and j terms 
scale similarly with £ provided that s/ + 3(sb + 1) = Sb~ 1. 

To summarize, in order for the scaled solutions, Eqs. (12 61 . 
to satisfy the Vlasov-Maxwell equations, it is necessary (and 
sufficient) for the similarity indices to satisfy 

s, = l, s p = s B +l, s f = -4-2s B . (32) 

In addition, we have shown that -1 < s B <Q and that 

LocD. (33) 

The relation between s p and sb implies that the Larmor ra- 
dius of the particles scales as L. The relation between s/ and 
sb implies that the energy density of accelerated particles in 
any momentum interval, Uh = XL / d 3 pcp f a , scales as the 
magnetic field energy density (when the momentum interval 
scales appropriately). Combining the above results, we find 
that (U h ) oc (U B ) oc Z> oc D 2s «. 

An illustration of the plasma configuration downstream of 
a hypothetical self-similar collisionless shock appears in Fig- 
ure Q. 




FIG. 1. — Schematic illustration of the model. The shaded vertical band 
at D = z = represents the shock. Shaded Circles designate regions where 
the magnetic field amplitude B exceeds some constant threshold. The charac- 
teristic scale L for variations in B increases with distance D from the shock. 
Curved lines illustrate the trajectories of two high-energy particles of dif- 
ferent momenta, p\ < P2, becoming magnetized (Rl < L) at two different 
scales. 
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3.5. Additional implications 

The scaling relations derived above may be used to con- 
strain the momentum dependence of the particle distribution 
functions. Using Lq = L(pq/ p) x l Sp in Eq. ( 1241 we obtain 



(/»>*,£)=(£) (fa 



PpoMf 



p/po— >oo 



(f a )(Wo,0) (*) 



l/.s,, 



«//»? 



(34) 



This implies that (f a ) (p,L) approaches a power-law distribu- 
tion at large p, (f a )(p,L) oc p'f/'p, provided that (f a )(p,L) 
approaches a finite limit for some p as L — ► 0. The condi- 
tion that (/ Q ) (p,L) — > Const as L — > is equivalent to the 
requirement that accelerated particles reach the shock front. 
This is physically reasonable since cosmic rays are expected 
to reach the upstream. Moreover, it is a necessary condi- 
tion for any acceleration process that involves the shock front 
(such as diffusive Fermi acceleration). 

A similar constraint may be obtained by taking the limit 
p/po ~ ¥ 0, which yields 



</ Q >(p,L)=(£ 



Sf/Sp 



(fa) 



VPo 



f/ s P 



p/po->0 



(f a )(ppo,oo)^Y \ (35) 



This implies that (f a ) (p,L) approaches a power-law distribu- 
tion at small p, (f a )(p,L) oc p s fl s '\ provided that (f a )(p,L) 
approaches a finite limit for some p as L — > oo. The condition 
that (f a ) (p,L) — > Const ^ as L — > oo is equivalent to the 
requirement that accelerated particles reach downstream in- 
finity. This is physically reasonable since cosmic rays are ex- 
pected to be advected with the flow to the downstream (e.g. in 
diffusive Fermi acceleration). Moreover, the presence of high 
energy electrons in the far downstream (~ 10 1( V J£ /) is required 
to account for afterglow observations. Thus, we conclude that 
the self-similar distribution function approaches a power-law 
dependence on p at both high and low momenta, 



dn a 
dp 



ocp 2 d 2 f>f a (pp) oc p 



-2/(^+1) 



(36) 



Let us consider next the magnetic field correlation function. 
Using L/Lo = Ax/Axq in Eq. (125 \ we obtain 



B, 7 (Ax,Af,L)=(^) 2 % 



2.S B 



Ax/ Axq— »oo 



(^) X(Axoi,0,0),(37) 



where i is a unit vector in the direction of Ax. This implies 
that Bjj(Ax, At,L) approaches a power-law behavior at large 
Ax, Bij(Ax,At,L) oc Ax 2 - !B , provided that By(Ax, At,L) ap- 
proaches a finite limit for some Axq as L, At — > 0. This is 
expected to be the case under the assumption that the accel- 
erated particles reach the shock front, since in this case the 
number of particles with Larmor radius greater than (some 
given) Ax approaches a constant as L — > 0, and such particles 
are expected to generate a correlation between the magnetic 
fields at points separated by Ax. 

Eq. I37i implies that the magnetic power spectrum has a 
component which scales as k~ 2 ~ 2sB . Assuming that the result 

Bij(Ax,0,L) w (Ax/ Axq) 2 "' Bij (Axox,Q,0) holds for Ax > 



77AX0, where r/ is some dimensi onles s constant, the power 
spectrum may be written as (see § I2.3> 



2 7 r^f ) (k J .,L)« f d^Bijix'^L) 



-ik_i_ -x^_ 



+ (Ax T 2s - 

/jjA.to 

= /i (kj_,7yAx ,L) 



dx 1 Xj^ ^ 



d6 Bjj(Ax x' ± ,0,0)e 



ik : -x' 



+ (Ax ) 



-2sbj^-2-2se 



du U 



1+2.S-B 



k ± r)Ax 



2tt 



d6 B ij (Ax &,0,<y)e~ ik± n , 



(38) 



where, u = k± x' ± . In the long wavelength limit, < k± <C 
(ryAxo) -1 , the exponent in f is approximately constant, such 
that f\ is nearly independent of k^. In addition, the lower 
limit in the second integral, k±i] Axq, may be taken to zero 
(recall that s B > -1). Thus, 2TrB (2D \k ± ,L) » f(i 1 Ax {h L) + 

k~ 2 ~ 2s " fxQux > Axo) where we assumed that the integral in the 
second term converges. As > — 1, for kj 
gleet /1 so that 



we can ne- 



Z?5f ) (k ± ,L)ocfe- 2 - 2sa . 



(39) 

The scaling of various other correlation functions, such 
as (fafp), may be deduced in the method shown above. 
For example, the arguments leading to the conclusion 
Bij(Ax,At,L) oc Ax 2 '* for Ax — > 00 also imply that 
B u (Ax,At,L) oc Af 2ss for At -» 00. 

4. THE "FLUID" COMPONENT, DOWNSTREAM 



As explained in 33.31 the assumption that L is the only rele- 
vant scale implies that the microscopic motion of the thermal 
particles (referred to as the "fluid" particles) is unimportant 
for the global solution. One option is that the current carried 
by the fluid is negligible, jp -C jh, in which case these parti- 
cles do not affect the EM fields. The second option is that the 
current carried by this population is considerable jp ~ /'/, ~ j 
(i.e. that it scales as the electric current carried by the accel- 
erated component), and that the thermal fluid is described by 
effective equations independent of the microscopic scales l sc i 
and Rijh = p,i ,c/eB. The first option is unlikely since, as we 
show in §4.21 the fluid is expected to be highly conducting 
and thus probably carries non-negligible currents. In this sec- 
tion we study the implications of the second option, jp ~ j'/,. 
We assume that the appropriate effective degrees of freedom 
are those of a single fluid, since the skin depth l si j and thermal 
Larmor radius Ruh have an important role already at the level 
of two fluid equations. 

In order to obtain equations which are independent of the 
downstream bulk velocity of the fluid, which can not scale 
with L, we consider the flow in the downstream frame (where 
the bulk velocity vanishes). Note, that since s t = 1, the self 
similarity relations, e.g. Eqs. i26\ . hold both in the shock 
frame and in the downstream frame (and in any other frame, 
as long as we set t = as the time when the shock was at 
z = 0). As the current carried by the fluid is non-negligible, 
j F ~ j, it must have a self similar structure. We thus assume 
that all the dynamically important fluid quantities have a self 
similar structure. In particular, the fluid velocity field u (in 
the downstream frame) scales as 



u'(x,f) = f"u(^). 



(40) 
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In 34.2l we show, that for L 3> Run the fluid can be consid- 
ered as infinitely conducting, i.e. 



E = --uxB. 

c 



(41) 



This, together with Maxwell's equation V x E = -c _1 <9,B, im- 
plies that L/T ~ u and thus s u = 1-^ = 0. In other words, 
the velocity fluctuations scale trivially. In order to relate this 
to the magnetic field we consider the force density applied 
to the fluid, ~ ~f 2 ri\m p u/T where T is the characteristic time 
scale for variations. The force has a contribution from the 
pressure, ~ P/L where P is the fluctuating part of the pres- 
sure, and from the EM fields, ~ c~ l jjB ~ Ub/L where Ub is 
the magnetic field energy density. If Ub is not negligible com- 
pared to P, requiring that the two force densities, ^Hn\m p u/T 
and c -1 jjB ~ Ub/L, scale similarly we find, using u oc L° and 
T oc L 1 , that Sb = 0. If Ub is negligible, P ^ Ub, the force 
density is dominated by P/L and sb can not be constrained 
by considering the force density scaling. However, P ^> Ub 
is unlikely. Our basic assumption that the accelerated compo- 
nent plays an important role in the solution suggests that the 
fluctuations in the fluid are driven by it. As the interaction 
between the accelerated component and the fluid is mediated 
through the EM fields this implies that they have a consid- 
erable contribution to the force density applied to the fluid. 
Moreover, the observational evidence that e# is not very small, 
eg ~ 10~ 2 - 10 _I , does not leave much room for P (which is 
smaller than Ub/^b) to be much larger than Ub- We conclude 
therefore that 

Su =s B = 0, (42) 



which, using Eq. (13 2\ . implies 



1. 



(43) 



Eq. 1421 and (I43i imply that U u oc Ub oc Uh, where U u , Ub and 
Uh are the energy densities of the fluid fluctuations, magnetic 
fields and accelerated particles, respectively. Implications re- 
garding the correlation funct ions of fluid quantities such as u 
can b e dra wn in analogy to § 13.51 

In 34. II below we study the scaling properties of the fluid 
equations of motion, under the simplifying assumptions of an 
ideal fluid and that corrections of order e# can be neglected. 
A closed set of equations is derived for the fluid component, 
and it is demonstrated that self-similar (sca lable ) solutions of 
the equations exist only for s u = sb = 0. In 34.2l we argue that 
under the assumption Ruh <C L, the fluid can be treated as 
infinitely conducting (as we show, the exact requirement is 
R L ,th/L<S:{v s T/L, L/w s T}). 

One consequence of the scaling indices derived above is 
a divergence of the energy associated with the self-similar 
components. A flat power-law energy-spectrum of acceler- 
ated particles, dn/dE oc E~ 2 [Eq. d36i with sb = 0], which 
is a natural consequence of the model and agrees with obser- 
vations and with linear DSA theory, implies that the energy 
density diverges at large momenta. Similarly, a flat magnetic 
power-spectrum of the form Bfj D \k± , L) oc k~ 2 [Eq. (I39t with 
Sb = 0] implies that the magnetic energy density diverges at 
large wavelengths (the same applies to the kinetic energy as- 
sociated with the velocity fluctuations). Poss ible remedies for 
the divergence problem are outlined in § 14.31 but resolving the 
problem is beyond the scope of this paper. 



4. 1 . Fluid equations 

For simplicity, we focus below on relativistic shocks. We 
comment at the end of this subsection on the application of 
the analysis to non-relativistic shocks. As mentioned above, 
we consider the flow in the downstream frame. 

The total electric current is j = j F + j h . The fluid current may 
be written as 



J,= £vxB-J 4 --LftE. 



(44) 



In order to guarantee conservation of charge we must include 



V ■E = 4ir(p h + p F ) 



(45) 



as an independent equation, where p/, is the electric charge 
density carried by the accelerated component, 



d P/q, 



(46) 



p>£ p eBL/c 



and pf is the electric charge density carried by the fluid. Con- 
servation of fluid energy and momentum is described by 



jFv, 



(47) 



where Tp V is the fluid energy momentum tensor, jp v is the 
fluid four current, and F^ v is the electromagnetic tensor. The 
fluid velocity u is defined so that T ' = in a local frame mov- 
ing with velocity u. These equations can be closed with the 
addition of an "equation of state," relating the different com- 
ponents of T^ v in the fluid rest frame. 

As the energy density of the fluid must scale as L°, the 
above equations are consistent only with a scaling in which 
the EM fields and Tp V both scale as L°, implying s u = sb = 0. 
However, since the equations may contain terms that can be 
neglected, and the scaling laws are likely to be relevant to fluc- 
tuations in (imposed on a constant background) rather 
than to its average value, different scaling exponents may also 
be allowed. In order to examine this in more detail, we make 
the simplifying assumption of an ideal fluid, i.e. T 1 ' oc <5 ,J in 
the fluid rest frame. This is the case if the distribution of the 
fluid particles is approximately isotropic in the fluid's local 
rest frame, which is reasonable as their Larmor radius is much 
smaller than the length scale for variations in the magnetic 
field. Under this assumption, equation J47I can be written as 

d, [ 7 2 (w + w)] + V • [7V0 + w)u] = E • j F + d t P (48) 



and 



7 2 (wo + w) E-j F + d,P 
z (o, + u- V)u+ ?-z u 



= -VP+p F E+-j F xB, (49) 
c 

where Pq + P and wq + w are the pressure and rest frame en- 
thalpy per unit volume, respectively, assumed to consist of 
a constant part, Pq ~ wq ~ r ^n\tn p c 2 ', and a fluctuating part, 
P, w. Here, 7 = (1 -m 2 /c 2 )~ 1//2 is the Lorentz factor associ- 
ated with the fluid's velocity fluctuations. The fluid equa- 
tions must be simplified by identifying terms that are neg- 
ligible, in order to determine the correct scaling laws. As 
explained above, we assume that the magnetic force density, 
c~ l \j F x B| ~ L~ l B 2 ~ L~ l eBjjniin p c 2 , makes a considerable 
contribution to the momentum balance [Eq. j49H. This im- 
plies P < eBl^n\m p c 2 ~ esPo- We assume that w/wq ~ P/Pq, 
which holds for example for any equation of state where the 
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enthalpy is a function of the pressure alone, in particular this is 
true for a relativistic equation of state, wq + w = 4(Pq + P). We 
find w < eB^ninipC 2 ~ e^wo- Comparing the first and the last 
terms in Eq. ( I49> . combined with the above result, L/T ~ m, 
we find that -f 2 u 2 /c 2 ~ and thus 7 w 1 and the fluid fluctua- 
tions are non-relativistic. We thus find that the thermal energy 
P and the kinetic energy U u ~ wqu 2 of the fluid fluctuations 
are related to the magnetic energy by P, U u < t/fl. Note that 
E j F ~ uB 2 /L [since £/p ~ (u/cB)(cB/L)], Hence, equa- 
tions J48I and i49i can be approximately written as 

w V-u = (50) 

and 

Wn 1 

-^(<9 f + u-V)u = -VP+p F E+-j F xB. (51) 

c l c 

In this case, we may eliminate the pressure from the equations 
by taking the curl of Eq. ( 15 1> . yielding 

^ {<9,V x u- V x [u x (V x u)]} = V x (p F E+ -j F x B). 

c z c 

(52) 

Equations @, (l28\-(E9t. iBTY d50l and (|52j, constitute a 
closed set of equations for the unknowns U,/ Q (for momen- 
tum p > ^ p eBL/c) and B. As can be readily seen, these equa- 
tions imply (and are consistent with) a scaling sb = s u = 0. 

The above analysis is also applicable to non-relativistic 
shocks. A few comments should, however, be added. In the 
non-relativistic case, the electric field is much weaker than 
the magnetic field, E — (u/c)B < (\ s /c)B < B. In addi- 
tion, it is expected that T ~ L/u -C L/c. Thus, the terms 
dtfa (for relativistic momentum) and q a E in Vlasov's equa- 
tion Eq. d28i . and the term c _1 9,E in Maxwell's equation 
V x B = 47rc~ 1 j + c~ 1 <9,E, may be neglected. Since the par- 
ticles of the thermal component are non-relativistic in this 
case, wo/c 2 may be replaced in equations (15 0> and d52l > by 
its non-relativistic limit, pF m (mass density). The assump- 
tion w/wo < P/Po, which for non-relativistic particles reads 
PFm/pFmO ^5 P / Po, is somewhat less justified but still reason- 
able. As mentioned above, these differences do not modify 
the scaling relations. 

4.2. Infinite fluid conductivity 

Here we argue that the assumption Ruh *C L implies that 
the fluid can be treated as infinitely conducting (as we show 
below the exact requirement is R^/L <C {\ S T/L, L/v s T}). 

The fluid velocity is approximately equal to the the veloc- 
ity of the proton component of the fluid. We now demon- 
strate that the electric and the magnetic forces acting on the 
proton component approximately cancel each other and thus 
lead to Eq. (I41> . Consider the proton fluid in a box of size 
< L during a time scale < T, in which the magnetic and 
electric fields are approximately constant in space and time. 
The proton momentum inside the box is ~ ~f s m p nul? . The 
time derivative of the momentum in the box originates from 
the electric force, f F ~ enEL? ~ en(L/cT)BL 3 , the magnetic 
force, /b ~ enu/cBL 3 , and the flow of momentum into the 
box, fp < L 2 n("f s m p v s )v s . The ratio of the momentum flow to 
the electric force is of the order of 

fp < lsm p w 2 cT ^av^T 

f E ~ eBL L ~ L L K ' 

Hence, the only force that can balance f E is f B . The forces 
are balanced only if E = — c _1 u x B. If the forces are not bal- 
anced, the proton fluid is accelerated. The time it takes the 



velocity u to reach the scale L/T (for which E ~ —c 'ux 
B), is very short, of the order of j s m p n(L/T)L 3 /enEL 3 ~ 
(RL,th/vsT)T <C T. Therefore, we may assume that the forces 
are approximately balanced at any given time, and the fluid is 
nearly infinitely conducting. 

4.3. Diverging Energy 

As explained in the beginning of this section, the vanishing 
of the magnetic scaling index, sb, implies logarithmic energy 
divergences. The divergence may be prevented if the scaling 
relations given in Eqs. ( I26t and d40l are approximate, rather 
than accurate. For example, high-order terms, oc l/L, l/L 2 ..., 
that were neglected in our analysis may introduce logarithmic 
corrections to the scaling relations. Alternatively, physical 
processes that were not included in the analysis, such as cool- 
ing, may limit the parameter range over which self-similarity 
is applicable, e.g. self-similarity may hold up to some cutoff 
momentum p max or cutoff distance D m „ r . 

Mod eling GRB ex ternal shocks ( Wa xmanlfi995l) and SNR 
shocks (Zhang 1993) implies that in both systems the maxi- 
mal energy of accelerated ions is limited by the age (size) of 
the system, whereas the maximal energy of the accelerated 
electrons is limited by energy losses. If the accelerated elec- 
trons play an important role in the evolution, self-similarity 
may break-down at some scale associated with electron cool- 
ing. Otherwise, the accelerated proton configuration may re- 
main self-similar beyond the electron-cooling scale. In this 
case, the accelerated protons have a time-dependent energy- 
cutoff and may (if their spectrum is sufficiently flat) modify 
the structure of the shock with time, rendering our assump- 
tion of a steady-state inaccurate (self-similarity may still ex- 
ist). We are not aware of observational constraints (on the 
shock thickness, say) that rule out this possibility. 

5. EXTENSIONS OF THE MODEL 

In the previous sections we have motivated the self- 
similarity assumption based on observational evidence that is 
relevant directly to the downstream of strong, non-magnetized 
shocks. It is possible, however, that in other related circum- 
stances, in which power-law distributions of energetic parti- 
cles are interacting with non-ma gnet ized MHD plasmas, sim- 
ilar scalable solutions arise. In § l5.1l we examine the possibil- 
ity that the hydr o-m agnetic structure of the upstream is also 
self-similar. In § !5.2l we discuss time-dependent, homogenous 
toy models, in which energetic particles with a power-law dis- 
tribution in momentum are added to a plasma. 

5.1. Upstream 

It is expected that the accelerated componen t precursor gen - 
erates turbulence in the upstream ( Blandford & Eichlei 1987). 
In the context of GRB's, it was recently claimed that the mag- 
netic field energy density in the upstream is larger than that 
typical to the typ ical ISM by at least 3 orders of magnitude 
dLi & Waxmatill2 00fi). which suggests that high energy parti- 
cles generate waves upstream of the shock. If this is indeed 
the case, it is likely that higher momentum particles are im- 
portant at larger distances from the shock, which suggests that 
the characteristic length scale for variations in the fields grows 
with the distance from the shock. Therefore, the EM struc- 
ture in the upstream may also have a self-similar character, 
although we do not have as strong observational evidence for 
self-similarity in the upstream as we have for such behavior 
in the downstream. 
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Consider the conditions many skin depths upstream of the 
shock. The equations governing the accelerated component's 
distribution function and the EM fields are the same as for the 
downstream. In particular we expect the relation between the 
scaling indices to be 



s p \ =s B \ + l\ 5/1 



-2s 



si ■ 



(54) 



whe re th e subscript 1 denotes the upstream. As explained 
in § 13.51 assuming that the accelerated component reaches 
the shock front, we expect the distribution function to have 
a power law momentum dependence with a power law index 
l p = Sf/sp. From the same considerations, we must have in 
the upstream l p = Sf\/s p \. We thus find that the scaling in- 
dices of the upstream and downstream are equal. In particular, 
for the expected value sb = 0, the energy in the magnetic field 
does not decay as we go farther into the upstream. This result, 
which cannot be true for arbitrary distances in the upstream, is 
related to the diverging energy in accelerated particles. Note 
also, that the value of B, which is constant according to our 
analysis in the upstream and far downstream regions, is ex- 
pected to change considerably in the shock vicinity (distance 
of order l s d where the self similarity breaks down), and thus is 
probably different in the upstream and downstream. 

Although the scaling relations are the same for the upstream 
and the downstream, the accelerated distribution function can 
have a qualitatively different spatial dependence. Note that 
for sufficiently high momentum the distribution is continuous 
across the shock. In particular, for a given momentum, the 
accelerated particle distribution function in the upstream is 
expected to d ecay exp onentially (if the acceleration process is 
DSA, see e.g. lKirketalJ2 000). while in the downstream it is 
expected to remain approximately constant. 

The analysis of the fluid component in the upstream is more 
complicated, as the magnitude of the (velocity, pressure) fluc- 
tuations and the amount of fluid heating are unknown. 

5.2. Homogenous configurations 

In this subsection we consider the development of magnetic 
fields in homogenous configurations as a consequence of the 
interaction of the thermal fluid wit h the acce lerated part icles 
(for r elated simulations see, e.g., iLuce k & Bell 2000j iBelll 
2004). Such configurations are much simpler to simulate and 
analyze than space-dependent configurations, and may give 
insights into dynamically important mechanisms. 

If spatial inhomogeneity has an important role in the accel- 
eration mechanism in collisionless shocks (as it does in the 
case of Fermi acceleration), it is likely that no acceleration 
takes place in homogenous configurations. In order to see the 
effects of an accelerated component on the dynamics, it would 
be interesting to perform homogenous simulations with an ac- 
celerated component (with a power-law distribution function) 
included in the initial conditions. 

In the following, we assume a homogenous configura- 
tion with a relativistic accelerated component with some 
anisotropic distribution function / Q (p) = g Q (/i)/?''', p > p m i„, 
included in the initial conditions, where /i = cos(p • z) (the z 
axis is chosen to represent the direction of the shock normal; 
the cylindrical symmetry around this axis is present in planar 
non-magnetized shocks and there is no reason not to include 
it in homogenous simulations). We assume that the fluid com- 
ponent may be described by single fluid equations for which 
there is no inherent physical length scale. We assume ini- 
tial neutrality, in the sense that any charge density (electric 



current) implied by f a is compensated by an opposite charge 
density (current) carried by the fluid. 

Quite generally, a non-isotropic distribution is unstable. 
The distribution function is expected to evolve from the initial 
unstable configuration to a stable (probably isotropic) con- 
figuration by the development of instabilities followed by 
dissipation. Particles with larger momentum will respond 
more slowly to the generated EM fields (which are limited 
in strength by the initial energy) and on longer length scales 
(perhaps of the order of their Larmor radius). Fluctuations of 
the EM fields on small scales are expected to decay with time 
( Gruzinov 200 1 a) while EM field generation at later times due 
to instabilities in larger momentum regimes is expected to oc- 
cur on longer length scales. It is thus likely that the magnetic 
field length scale grows with time. Self similarity is to be 
expected once the length scale of the magnetic field is much 
larger than the Larmor radius of particles having Pm„. 

The time development may be affected by the choice of g a . 
In particular, the value of the bulk electric current carried by 

the initial population oc J [Agpip)- geip)\ and bulk velocity 

oc J_j lAgpi^)+ge(lA] may affect the nature of the instabilities 
involved. As we show in this subsection, the self similarity as- 
sumption (as long as it is valid) allows determination of some 
physically interesting quantities without dependence on the 
precise instability mechanism. 

In this case, we cannot assume a steady-state. However, 
there is a full (three dimensional) translational symmetry (ho- 
mogeneity). We consider averaged quantities at fixed times. 
In particular, Eq. (1151 and Eq. (1161 are replaced by: 

(/«)(P,0= lira -j—r / a! 3 x7 a (x+x',p,f) (55) 

r->oo J x , <r 



and 



fi, 7 (Ax,Af,r) = 

lim -j-—; I c/Vb^x+x' ,t)Bj(x + Ax+x' ,t + At), 

(56) 

where these functions do not depend on x due to the homo- 
geneity {B defined similarly). 

We assume that a scaling property, following Eqs. 1231 - 
d25l . is obtained at late times, i.e. 



(57) 



B(L) = ( — ) B(Lo), 



(fa)(V,L)-- 



Lo) ^ I. : 



and 



Bij(Ax,At,L) = 



L\ 2SB 



At 



Ax 



L/L Q : (L/L yi 



(58) 



(59) 



with L oc f'/ s ', Toe t (s, > 0). 

An important difference with respect to the shock case 
arises from the absence of a finite bulk velocity, and therefore 
the assumption that time scales as distance no longer applies 
in general. We study cases in which the magnetic fields are 
stronger than electric fields and thus restrict s, > 1 . A value 
of s, > 1 is consistent with Vlasov's equation, Eq. d28t . if we 
neglect the time derivative term and the electric field term, 
which is justified for s, > 1 . 
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The distribution function of the accelerated component is 
by definition finite at t — > (also at L — ► 0). In analogy to Eq. 
J34l . we therefore have (/ Q ) (p,f) oc for p — > 00. On the 
other hand, (/ a ) (p,f = 0) oc p which implies that S//s p = 
l p (this also shows that a finite initial distribution must be a 
power-law in order to achieve self-similarity). From Vlasov's 
equation (and the expression for the electric current) we find, 
as for the shock case, that 



s p = s B +l, s f = -2s B -4. 



(60) 



Here it is useful to solve for the scaling indices in terms of the 
initial power-law index l p , 



lp + 4 

'lp + 2' 



l p +2 



2lp_ 
'lp + 2' 



(61) 



Assuming t hat th e fluid can be considered as infinitely con- 
ducting (see § 14.21 . we have u ~ L/T so s u = 1—5,. Assuming 
further that the magnetic force is not negligible in the fluid 
momentum equation (15 1 1 we find s u = s B (this is reasonable 
as the energy source is the energetic component, and energy 
is transferred to the fluid through the magnetic fields). To- 
gether with Eq. J61 1 we have 



s, = l-s B = 



20, + 3) 
lp + 2 



(62) 



In particular this implies B oc L SB oc t SB ^ s ' = r^V^+S). 

6. SOME GENERAL IMPLICATIONS 

If self-similarity holds, it has important implications for any 
model of particle acceleration and/or field generation. In § 16.11 
we show that a previously suggested model, that d escribes the 
flow in terms of merging current filaments ( Medvedev et al. 
2005) in a self similar manner, must be generalized in order to 
be a self-consistent model. We show that after the appropriate 
gen erali zation the model follows the scaling relations derived 
in § 15.21 and that the gen eraliz ation significantly modifies the 
model's predictions. In § 16. 21 we discuss the relevance of our 
analysis to DSA. 

6.1. Current merger model 

Med vedev et al J (12005 ) suggest a model of coalescing elec- 
tric current filaments in a quasi two-dimensional configura- 
tion for describing the magnetic field evolution in the down- 
stream of collisionless shocks. In this model there is no dis- 
tinction between accelerated and thermal particles, so it does 
not hold information regarding the energy dependence of the 
particle distribution function. As the model is essentially ho- 
mogenous and self-simi lar, it is interesting to compare it to 
the scalings derived in § 15.21 

This model assumes that electromagnetic instabilities in the 
shock lead to the formation of current filaments, which may 
be approximated as infinite in length and having equal electric 
currents oriented parallel to the flow, half of the currents posi- 
tive (oriented along the flow) and half negative. The filaments 
are assumed to be distributed randomly in space. This sim- 
ple model assumes that current filament coalescence evolves 
in discrete steps. In each step the filaments carrying similarly 
oriented currents, are divided into neighboring pairs, which 
attract each other and merge. Two filaments with diameter D, 
inertial mass per unit length fi, and current /, unite to form a 
new filament of diameter D' = \/2D (conserving area), mass 
per unit length // = 2fi, and current /' = 21. In each step, the 



typical distance d between neighbor filaments with similarly 
oriented currents, grows by a factor of v2, d' = \[2d, as a 
consequence of the reduced filament number density. It is as- 
sumed that d ~ 2D which implies that the filaments roughly 
fill the space with their area. Initially, all the filaments are 
identical (except for the sign of the current), and are hence 
identical at any given time. The coalescence time, here de- 
noted by T coa , is estimated by analyzing two isolated, parallel, 
infinite current filaments attracting each other. 

The fact that D and d change in each merger by the same 
factor, implies that the configuration following each merger is 
identical to that before the merger, with re-scaled parameters. 
In particular, the filament coalescence time, t coo oc fj}' 2 cD/I, 
and the maximal velocity during coalescence, v oc //(c/x 1 / 2 ), 
scale as t' coo = r coa and v' = \f2v. The fact that the filament ve- 
locity grows with time, led the authors to the conclusion, that 
velocities of the order of the speed of light will be reached. 
Once the filaments move with v ~ c, the coalescence tempo- 
ral behavior changes considerably. 

This suggested merger model is problematic, as it implies 
that the magnetic field grows by a factor of V2 in each step. 
For example, the magnetic field produced by a single filament 
at its edge, B = 4/ /cD, scales as B' = \/2B. More generally, 
the coalescence suggested is self-similar with j'(x) = j(x/\/2), 
implying B'(x) = %/2B(x/V2). To see this, note that the cur- 
rent density carried by each filament, j ~ I/D 2 , does not 
change in the coalescence, while the filament diameter and 
distance between filaments grow by y2. The fact that the 
magnetic field grows by a constant factor in each step, im- 
plies that after a few steps the magnetic energy density grows 
beyond equipartition, which is impossible. This can be cor- 
rected by changing appropriately the coalescence conditions. 

Since the physical process of the merger is complicated, 
there is no simple way to determine the post-merger cur- 
rent directly. In general, one should therefore set the cur- 
rent of the merged filament to {\frfl where £ is a parame- 
ter of the model (the mass per unit length is set to 2fi since 
mass is conserved). In this case, the coalescence is self- 
similar with j'(x) = (v / 2)''~ 2 j(x/v / 2), implying that B'(x) = 
(v / 2)^~ 1 B(x/v / 2) and the magnetic field amplitude does not 
grow for ( < I . The scaling of the filament coalescence 
time and maximal velocity changes to t' coci = (v2) 2 ~^T coa and 
v' = (V2) c_1 v. 

The (necessary) generalization of the current merger model 
strongly influences the conclusions that can be drawn based 
on this model. The growth in length scale is only determined 
up to a free parameter £. In order for the magnetic field not to 
diverge with time, we must have £ < 1 (if the magnetic field 
amplitude is postulated to be constant in time, ( = 1). This 
implies that the velocities do not grow (remain constant for 
£ = 1) and therefore do not reach the speed of light. 

The (corrected) scaling s ar e compatible with the scalings 
for the fluid derived in § 15.21 with Sb = (— l,s t = 2— £,s u = 
£ — 1 . Note, that the self-similar analysis allows us to reach 
most of the physically interesting conclusions without making 
oversimplifying and model-specific assumptions, such as the 
calculation of the length of the time step using two isolated 
filaments in the current merger model. 

It has been claimed (Milosavli evic & NakaJ 12005 al) that 
current filaments are unstable to kink-like modes and that 
the two dimensional configuration is therefore disrupted. We 
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note that even if the instability is efficient, this does not rule 
out the model. Since the instabilities are derived within the 
framework of MHD, for which the above scalings hold, the 
perturbation's growth rate scales as all other time scales, i.e. 
r' im = {\f2) 1 ~ C 'Ti ns where r,„ s is the inverse growth rate. If ini- 
tially Tins > T coa (which is not ruled out), the filaments merge 
before they are destroyed by the instabilities, and this holds, 
i.e. > T coa, for all subsequent mergers. 

6.2. Diffusive shock acceleration 

Diffusive shock acceleration (DSA) is the mechanism be- 
lieved to be responsible for the production of non-thermal 
populations of charged particles in collisionless shocks. In 
this first-order Fermi acceleration process, particles gain en- 
ergy by repeatedly bouncing between the converging up- 
stream and downstream flows. In the lack of a self-consistent 
theory for the interaction between electromagnetic waves and 
accelerated particles in collisionless shocks, most progress 
was made under the "test-particle" approximation. In this 
approach, the effects of the waves are modelled by some 
particle-scattering Ansatz, while the influence of the parti- 
cles on the waves and on the shock structure are neglected 
(for reviews, see iDrurvl 119831 iBlandford & Eichlerl 11987^ 
iMalkov & Drur\l2001l) . 

In the case of non-relativistic shocks, DSA has been 
successful in reproducing the power-law spectra observed 
in strong shocks, under very general assumpt i ons regarding 
the scatte ring mechanism ([Kr vmskiil f 1 9771 lAxford et al.l 
fl977t lBelH978t [Blandford & Osrtkeilll978l) . The analysis 
of relativistic shocks is more complicated, mainly because 
the non-thermal particle distribution is not isotropic. The 
particle spectrum was calculated in relativistic shocks under 
various assumptions regarding the scatter ing mechanism , 
using mon t e-carlo simulations [e.g. lEUison et alJ J l99Q): 



from the shock, l p < -4 would imply that the currents carried 
by the accelerated component j h decrease faster than the to- 
tal current j ~ V x B, and thus the effect of the accelerated 
particles on the magnetic field is negligible at large distances 



Ostrowski ( 199 11)1, and by numerical jK irk & Schneider 
1987 tlGallant & Achterberdl999tlvietril2003UBlasi & Vietril 
20051) or analytical llKeshet & Waxmanl 120051) study of the 
transport equations. In general, a power-law distribution 
function of accelerated particles is found, with indices that 
depend on details of the model and usually satisfy l p < -4. In 
the case of ultra-relativistic shocks with isotropic, small-angle 
scattering, numerica l studies have converged on a spectral 
index l„ ~ -4.22 llBednarz & Ostrowskil 119981 iKirk et all 
1200(1 lAchterberg et al.1 120011) . in agreement with GRB 
afterglow observations and with the analytic result l p = 38/9 
l lKeshet & W axman 2005) . However, this value was demon- 
strate d numerically ( Ballard & Heavens 1992; Ostrowski 
19931 lEUison & Double! 120021 120041 iMeli & O uenbv 
2003b a Bed narzl 120041 iNiemiec & Ostrowskil [ 20041: 
Lemoine & Pelletierl 120031 iLemoine & Revenul 120061) and 
analytically (Keshet & Waxman 2005) to be sensitive to the 
scattering mechanism, which is poorly constrained. 

In our analysis, the fluctuations in the fluid provide a natural 
scattering mechanism for the accelerated particles. Since the 
analysis reflects near equipartition between fluid fluctuations, 
accelerated particles and magnetic fields, it does not naturally 
evoke (but does not rule out) a test-particle approach. A spec- 
tral index l p jl 4 can be reconciled with our analysis only if s b 
we relax some of our assumptions. For example, if we as- 
sume that the magnetic field energy is constant with distance 



(in this case, the test particle assumption is self-consistent). 
If, on the other hand, we assume that the magnetic field de- 
creases with distance as sb = —(l p + 4)/ (l p + 2) (see § 13. 5> . a 
value l p < -4 would contradict our assumptions regarding the 
self similarity of the fluid component (for example, the mag- 
netic force would become negligible in the momentum con- 
servation equation, or the fluid would no longer be highly 
conductive and its currents would become negligible). 

The assumption of self-similarity places constraints on 
some properties of DSA. As an illustration, consider the case 
of small-angle scattering, parameterized by a propagation- 
angle diffusion coefficien t D,,.,,.. The stationary transport 
equation can be written as (Kirk & Schneid erll987l) : 

ji(Vi/c+fj,i)d Zt (f) (pi,fjLi,z) = d m [D^(jn,pi,z)d m (/,■)] 

(63) 

where subscript i = 1,2 denotes upstream or down- 
stream parameters, respectively, 7/ = (1 - v?/c 2 ) -1 / 2 , and 
[i = cos(p,- • z). The momentum p ; (and therefor also 
/i) is measured in the fluid rest frame, whereas z is 
measured in the shock frame. With boundary condi- 
tions {fi)(pi,Hi,-oo) = 0, {f) (p h m, 00) = (/,■) ^(pi) and 
(/i)(Pi»Mi)0) = (/2)(p2,M2,0), where (pufJ-i) and (P2,P>2) 
are related through an appropriate Loren z boost. The solu tion 
is unique up to the normalization of (/,) (Kirk et al. 2000). 

The small-angle scattering assumption requires that the 
Larmor radius of the particles Ri is much larger than the 
magnetic field length-scale L. In our framework, for any 
given momentum p this is true up to a limited distance from 
the shock (where L < pc/eB). Consistency of DSA theory 
thus requires that (f)(pi,Hi,z) converge to its limiting value 
(fdoo within this range. 

In order to comply with the self similar scaling of Eq. i24l . 
we must require that 



(fi)(PhHi,z)-- 



Z() 



(fi) 



(z/zo)' 5 " 



,P>i,Z0 



(64) 



This scaling can be reconciled with the transport equation 
[Eq. (1631 1 if and only if the diffusion function scales as 



(z/zo) s " ' 



M;,Zo 



(65) 



The boundary conditions are invariant to the above scaling (at 
the shock front p\ oc pi). Under this scaling, (fi) (Pi, p,z) oc 
p 1 '' at z = ±00, where l p = sj/s p . The values of Sf and l p are 
determined in a non-trivial way by the function (/?,-, /Lt,-,z). 
If s p = 1, as is expected when the plasma is highly conductive, 
we may write Eq. d65i as 

Df_ lfl (p i ,Hi,z) = z~ l D^(Hi>Pi/z). (66) 

As an example of the above scaling, consider a highly rel- 
ativistic particle scattered by weak magnetic fluctuations of 
correlation length L <C Rl- The particle trajectory may then 
be described as a random walk in fi, and the resulting dif- 
fusion function is D^^ip, p,,z) = Di^LRJ 2 . For L oc z (and 
0), we obtain D w (p,^,z) oc z~ l {p/z)~ 2 . 



7. DISCUSSION 



We have studied the consequences of the assumption that 
the downstream flow of non-magnetized, collisionless shocks 
is self-similar. This assumption is motivated by the existence 
of a strong magnetic field many skin depths downstream of the 
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shock front, as inferred from observations of GRB afterglows 
and young SNR's. This suggests that the correlation length 
of magnetic field fluctuations, L, diverges with the distance D 
from the shock front. Although our analysis was motivated 
by evidence for self-similarity in GRB afterglows, and possi- 
bly also in SNRs, it may apply to other systems with similar 
characteristics, such as the shocks involved in the large scale 
structure of the Universe. 

As the electromagnetic fields and particle distributions at 
any given point fluctuate with time, a stationary shock struc- 
ture may be described by the averages (over planes parallel 
to the shock front) and correlation functions of the fluctuating 
quant ities, which depend only on the distance from the shock 
(see § I2.3> . The self-similarity assumption implies that L oc D 
and that the averages and correlation functions at different dis- 
tances D from the shock, corresponding to different values of 
L, are rela ted to each other by simple scaling transformations 
(see §E3, e.g. 



u:.)(v-L)= [ ^) f (f a ) 



B u (Ax,At,L) = 



2.v B 



Bi 



At 



Ax 



L/Lo' (L/L y 



-,L 



where the scaling exponent of the magnetic field must satisfy 
-1 < s B < 0. 

A schematic illustration of a self similar downstream config- 
uration is presented in figure [2 

We have argued (see § 13.31 that the similarity assumption 
suggests that the plasma may be approximately described as 
a combination of two self-similar components: a kinetic com- 
ponent of energetic particles, and an MHD-like component 
representing "thermal" particles. We have argued that the en- 
ergetic particles are likely to carry a significant fraction of the 
current, and derived (using the Maxwell- Vlasov equations) 
the scaling of the characteristic time for variations in the phys- 
ical quantities, the scaling of the characteristic particle mo- 
mentum and the scaling of the particle distribution function 
normalization: 



s, = 1 , 



■s B + l, s f = -4-2s B . 



These relations imply that the characteristic Larmor radius of 
energetic particles scales as L, and that the energy density of 
energetic particles in any momentum interval, with the inter- 
val scaling as Z>, scales as the magnet ic fi eld energy den- 
sity oc L 2s ". We have then shown (see § I3.5> that (under the 
assumption that accelerated particles reach the shock front 
and/or are advected to the downstream) the spectrum of ac- 
celerated particles is 

dn/dEocE- 2 ^ SB+l \ 

and that the scaling of the magnetic correlation function (for 

Ax — > oo) is 

(Bi(x)Bj(x+ Ax)) cx Ax 2sa . 

Similar conclusions can be drawn regarding various other cor- 
relation functions. 

The "thermal" particles were discussed in § 0] In § I4.2l we 
have argued that the thermal component may be considered 



as an infinitely conducting fluid. We have shown that in this 
case sb = and the scalings are completely determined, e.g. 
dn I dE oc E~ 2 and B oc D°, with possible logarithmic correc- 
tions. We have derived in § !4.1l a closed set of equations for 
the fluid component, under the simplifying assumptions of an 
ideal fluid and neglecting corrections of order Eb- 

The self similarity assumption and its implications do not 
hold for arbitrarily large distances and high particle momenta 
where new physical processes that were not taken into account 
in our analysis come in to play. For example, our assump- 
tion of an infinite planar shock is invalid at distance scales of 
the order of the blast wave radius, and the use of the Vlasov 
equation is invalid for high momenta where radiative effects 
cannot be ignored. Upper cutoffs to the distance scale and the 
momentum range described by the self-similar solution, or 
logarithmic corrections to this solution, are possible re medie s 
of the energy divergence when sb = 0, as discussed in § 14.31 

If self-similarity holds, it has important implications for any 
model of particle acceleration and/or field generation. In § 16.21 
we have shown that the velocity-angle diffusion coefficient for 
small-angle scattering in diffusive shock acceleration models 
must satisfy D^^p^D) = D~ l D^^(\)/D) (where p is the parti- 
cle mo mentu m). 
In § 16.11 we have discussed 



the 



model suggested by 
iMedvedev et alJ J2005I) for the generation of a large scale 
magnetic field through hierarchical merger of current- 



filaments. We have shown that in order to avoid a diverging 
magnetic field, the model must be generalized by allowing 
a more general scaling of the electric current of merged fil- 
aments. The gen erali zed model follows the scaling laws we 
have derived in § 15.21 This implies that the self-similar anal- 
ysis allows us to reach most of the physically interesting con- 
clusions without making oversimplifying and model-specific 
assumptions (such as those related to the calculation of the 
merger time). The predictions of the generalized model dif- 
fer substantially from those of the original model. It predicts, 
e.g., a scale independent merging velocity, rather than an in- 
creasing velocity approaching the speed of light (this is valid 
for a non-decaying magnetic field; for a decaying field, the 
velocity decreases with scale). Finally, we have shown that 
the instability of current fi laments (which was pointed out by 
Milosavlievic & Nakar 2005a) does not necessarily imply that 
the cu rrent merger model is not viable. 

In § !5.2l we have pointed out that the self-similarity assump- 
tions may be tested through their predictions for the evolution 
of homogenous (time-dependent) plasmas, which may be ac- 
cessible to direct numerical simulations. An inclusion (at the 
initial conditions) of an anisotropic, power-law spectrum of 
high energy particles, dn/dE oc E 2+l '', in homogenous simu- 
lations may lead to a self-similar evolution in time, described 
by the scaling exponents given in Eq. d61t . Assuming that the 
fluid can be considered as infinitely conductive, we predict 
that the magnetic field evolution follows B oc f-( , p +4 )/ 2 ('fj+ 3 ). 

Our self-similar model predicts the scaling of all physical 
quantities related to the accelerated particles, electromagnetic 
fields and fluid fluctuations. In particular, the particle spec- 
trum is related to the magnetic field scaling: An accelerated 
particle distribution with a spectral index l p < -4 indicates a 
decay of the magnetic field amplitude with distance from the 
shock according to £r('p +4 )A'/> +2 ). Such a decay of the mag- 
netic field may be detectable, for example through the spatial 
dependence of synchrotron emission from non-thermal elec- 
trons gyrating in the downstream magnetic fields. 
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Finally, we should emphasize two major open questions re- divergence is identified in the energy of magnetic fields and 
lated to our analysis, which need to be resolved. First, the va- fluid fluctuations, 
lidity of the description of the "thermal" component in terms 
of single-fluid MHD equations requires verification. Second, 
the diverging energy in accelerated particles that results from 

the distribution n(E) oc E~ 2 , indicates that our assumptions This work was supported in part by ISF and AEC grants, 

cannot be valid for arbitrarily high momen tum (see discus- We thank A. Gruzinov, A. Loeb & A. Sagiv for helpful dis- 
sion and possible remedies at the end of § 13.51 . A similar cussions. 

APPENDIX 

COLLISIONLESS SHOCKS IN 3D PARTICLE-IN-CELL SIMULATIONS 

In recent years, particle in cell (PIC) simulations have provided important preliminary clues to the nature of collisionl ess shocks. 
Fully three-dimensional PIC simulations have began to explore as pects of such shocks in non-homogen eous ( Nishikawaetajj 
l2003l fFrederiks en et al.l2004lJSpitkovskvl2005l) and homogeneous l lSilva et al.l2003tlJaroschek et alJ2004llRomanov et alJ2004) 
flows. The main conclusions relevant for collisionless shocks are: 

• Confirmation that shocks in pair-plasma (with no highly accelerated particles) are mediated by electromagnetic, two-stream 
(Weibe l-like) in stabilities if the upstream is sufficiently non-magnetized (magnetic energy less than ~ 1 % of the kinetic 
energy; Spitkovskv 2005) and cold. In interpenetrating shells, this instability produces filaments in density (if T > 5) and 
in electric current (Jaroscheket al. 2005), initially with <skin depth thickne ss, and transve rse (perpendicular to the flow) 
magnetic fields. Electrons in ion-electron plasma behave similarly dFrederiksen et alJ2004l) . 

• Prelim inary evidence for the existence of relativistic, collisionless shocks in a non-magnetized pair-plasma. Spitkovskv 
(2005) resolves a shock-like behavior in a 3D simulation of interpenetrating pair-plasma shells, each with Lorentz factor 
r = 15. A density jump in approximate agreement with the Rankine-Hugoniot adiabat extends over ~ 70/«/, and coincides 
with a region of near-equipartition magnetic fields. Here, l sl j = c/ui p and ui p = (Aitne 1 /Tm) l l 2 are the skin-depth and plasma 
frequency, respectively, of simulated particles of mass m, number density n and Lorentz factor T. 

• Possible indications for self-similarity. In interpenetrating plasma shells, the filaments merge after saturation and grow in a 
hierarchical process, maintaining a configuration of similar featu res on gradually la r ger scales, with m agnetic fields decay- 
ing slowly and possibly saturating at e B < 1% [in pair plasmas: ISilva et alJ (120031) : iMedvedev et a Frederiksen 
et al. (2004; Figure 3) provide evidence that the power-spectrum of transverse magnetic fields is roughly a power-law, 
Pk oc k~ a with a ~ 2.3 - 3.0, extending to gradually larger scales. 

Present 3D PIC simulations are limited to small simulation-box volumes, V < i0 6 l^ d , and to short simulation durations, T < 
lOPwZ 1 . Here, u> p and l s d corresponds to the species with fastest-response; simulating an ion-electron plasma implies much 
smaller volumes and durations when written in terms of the ion parameters. In order to resolve any effects associated with the 
ions, an effective, small proton to electron mass ratio, m p /m e < 20 (with present computational resources) must be used, and the 
preliminary results thus obtained are not easily extrapolated to more realistic mass ratios. These constraints limit the relevance 
of any conclusion drawn from such simulations, and may even undermine their reliability: 

• The length of the simulation box parallel to the flow is too small in most cases to resolve a shock , even in pair-plasma 
[with the possible exception of Spitkovskv 2005]. The configuration found in most simulations thus represents merely a 
transient stage in the formation of a shock, which may bare little relevance to its steady-state. Moreover, the short length 
of the box suppresses longitudinal modes of long wavelengths, thus distorting the plasma evolution by effectively reducing 
it to 2D. For example, an expo nential decay of the magnetic energy (after growth saturation) abruptly stops in simulations 
(ISilva et al!2003ll.Taroschek et al.l2 005). possibly due to this effect. 

• The small size of the simulation box (< tens of skin-depths) perpendicular to the flow places an artificial cutoff on transverse 
modes with long wavelengths. Simulations indicate that the transverse scale o f the most energetic modes grow s rapidly, 
reaching the box size cutoff at early stages of the simulation (Silva et al. 2003, Figure 3; Frederiksen et al1 l2004l Figures 
3-4). After this occurs, the simulations are strongly affected by the boundary conditions, and the simulated evolution is 
probably highly distorted. 

• The Larmor radius of particles of Lorentz factor 7 is Ri = l s d[tB(l 2 - l)/2r(r- 1)]~ I//2 , where V is the average Lorentz 
factor. A simulation box tens of skin-depths long is therefor not sufficiently large to simulate particle magnetization when 
the ratio between magnetic and thermal energy densities is eg < 1%. Moreover, the box is much too small to resolve 
acceleration of particles to high energies, so any Fermi-lik e acceleration process is su ppressed. Electrons are observed to 
be accelerated in ion-electron simulations (Nishikawa et al.l2003UHededal et al.l20 04). but to energies smaller than Tm p c 2 . 
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